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A hydrodynamic-type, macroscopic theory was set up recently to simultaneously account for 
dissipation and dispersion of electromagnetic field, in nonstationary condensed systems of nonlinear 
constitutive relations Since it was published in the letter format, some algebra and the more 
subtle reasonings had to be left out. Two of the missing parts are presented in this paper: How 
algebraically the new results reduce to the known old ones; and more thoughts on the range of 
validity of the new theory, especially concerning the treatment of dissipation. 
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I. INTRODUCTION 

The macroscopic Maxweh equations, given in terms 
of E, D, H, and B, need constitutive relations linking 
them to be closed. The form of these relations depend 
crucially on two physical parameters, the frequency and 
the field strength. Weak fields are necessary for the valid- 
ity of the linear response theory, E ~ D,H ^ B; while 
the hydrodynamic constitutive relations presuppose 
small frequencies. 

While the limits on the constitutive relations are well 
known, the equally important closure problem one level 
higher has not entered the general awareness: A elec- 
tromagnetic theory is only complete if it also consid- 
ers the feedback, the force on the volume elements of 
condensed matter, exerted by the electromagnetic field. 
In the microscopic electrodynamics, this is simply the 
Lorcntz force - the Maxwell equations account for the 
field produced by the sources, while the Lorentz force (in 
conjunction with the Newton equation) describes how the 
field changes the position and the motion of the sources. 
In the macroscopic theory, two quantities are necessary 
to close the theory: (i) The additional energy due to the 
presence of the electromagnetic field, and (ii) the flux 
of the conserved, total momentum density, ie the total 
stress tensor including both the material and the field 
contribution. The hydrodynamic theory provides unam- 
biguous expressions for both, and is therefore closed and 
complete. Circumstances are less fortunate for the lin- 
ear response theory, as these expressions are only known 
with a string of additional restrictions. 

Assuming transparency (ie lack of dissipation) , quasi- 
monochromatic external field and stationarity (ie iden- 
tically vanishing velocity field of the condensed system) 
— all in addition to the linearity of the constitutive rela- 
tions — Brillouin obtained the field energy in 1921, while 
Pitaevskii, forty years later, arrive at the attendant ex- 
pression for the total stress tensor, cf §80, 81 of the classic 



book of Landau and Lifshitz |^ , and the comprehensive 
and informative revew article by Kentwell and Jones [Q. 

If we draw a diagram of field strength versus frequency 
cj, with the field strength pointing to the right, and lo 
upward, a vertical stripe along the w-axis represents the 
range of validity for the linear response theory, while the 
hydrodynamic theory reigns within a horizontal stripe 
along the field axis. The expressions of Brillouin and 
Pitaevskii are valid in isolated patches in the vertical 
stripe, wherever dissipation is negligible Q. 

The parameter space beyond the above two perpen- 
dicular stripes needs a theory that can simultaneously 
account for dissipation, dispersion, nonlinear constitu- 
tive relations and finite velocities. Although one might 
expect principal difficulties in setting up such a theory, 
as neither of the two parameters, field and frequency, re- 
mains small, we are (up to and maybe slightly beyond 
the optical frequencies « lO^^Hz) still in the realm of 
macroscopic physics, as the electromagnetic wavelength 
remains large compared to the atomic graininess. And 
when asking questions such as what is the force on a vol- 
ume element exerted by a strong laser beam, if we confine 
our curiosity to the averaged force — with a temporal 
resolutions larger than the time needed to establish local 
equilibrium — a simple, universal and hydrodynamic- 
type theory is still possible. Such a theory was derived 
recently, which includes the dynamics of polarization, but 
(in this first step) neglects magnetization [Q . 

In this paper, we show explicitly that the new and 
nonlinear expressions for the energy and the stress ten- 
sor indeed reduce to the known ones, of Brillouin and 
Pitaevskii, in the specified limit. Because the non- 
linear theory is the result of a qualitatively different, 
hydrodynamic-type approach, and not simply a "gen- 
eralization" of the linear theory, this outcome is by no 
mean obvious and assured. Besides, the associated alge- 
bra is fairly involved and needs to be presented. Once 
accomplished, this provides two additional bonuses: It 
strengthens our trust in the new theory and provides a 
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transparent interpretation for the old and classic results 
which, having relied heavily on algebra, are somewhat 
lacking in appropriate physical pictures. 

The hydrodynamic theory of dispersion will be pre- 
sented in section Q to render the present manuscript 
self-contained. Because we shall only list the relevant 
formulas and abstain from repeating all the reasoning 
and arguments that lead to the new theory, the reader is 
advised to also read However, the important ques- 
tion on the range of validity of the new theory in dis- 
cussed here in greater details than in [Q, at the begin- 
ning of the next section. Section [I] incorporates the 
specified approximation and deduce four results. They 
are compared to the energy density by Brillouin, and to 
three formulas by Pitaevskii's: the total stress tensor, 
the "nonmagnetic" magnetization, and the time depen- 
dent permittivity. Section ^ ends with a brief summary. 



II. THE HYDRODYNAMIC THEORY OF 
DISPERSION 

In this section, we shall first discuss in some details 
the range of validity of the new theory, then present its 
complete set of equations, and specify the theory by an 
expansion of the thermodynamic energy to third order in 
the field variables. 



A. The Range of Validity 

The proper hydrodynamic theory of electromagnetism 
accounts for the macroscopic dynamics of continuous 
media in the low frequency limit, for a system that is 
charged or exposed to external fields. Local thermody- 
namic equilibrium holds, and the set of hydrodynamic 
variables is identical to that of the thermodynamic vari- 
ables. The equations of motion are conservation laws 
and the Maxwell equations, including irreversible terms 
accounting for dissipation. 

At higher frequencies, microscopic variables deviate 
more and more from equilibrium, becoming independent, 
and finally ballistic. Denoting the time tioc needed to 
establish local equilibrium, this starts to happen when 
UTioc is no longer small. To account for this circum- 
stance, we usually have to abandon the hydrodynamic 
theory and embrace the Boltzmann theory which, despite 
its undeniable usefulness, is both a far more complex and 
a less general theory - it considers the vast number of mi- 
croscopic degrees of freedom explicitly, and it is confined 
to dilute systems. The question therefore is whether an 
appropriately generalized hydrodynamic theory can be 
made to account for some of the more interesting aspects 
at higher frequencies, and save us from the complexities 
of the Boltzmann theory. 

Let us concentrate on one specific microscopic variable, 
the polarization P. Actually, as far as its spatial extent is 



concerned, it is a macroscopic rather then a microscopic 
variable, but it is certainly dependent in the hydrody- 
namic limit, as long as ujTp <^ 1, where rp is P's longest 
time scale. (All this can also be said of the magnetiza- 
tion, which we however shall not consider here.) 

In a dielectric medium, P has many characteristic 
times, given by the resonance frequencies and their 
widths. If they are well separated, then the equation of 
motion, close to one resonance and in the simplest case 
considered below is 



F/luI - tP + P = xD, 



(1) 



cf Eq(|lJ) of the next section. The given resonance 
may be overdamped (r ^ l/'^p) or sharply resonating 
(r <C l/(jjp), and the characteristic time rp, after which 
the polarization is no longer independent, is respectively 
r and 2/{ujpT), while l/wp is typical for the time scale of 
P's motion in the resonating case. (Note that going to a 
different resonance, rp will change, it therefore depends 
on the frequency of the external field.) 

If the polarization P is a specially slow variable, rp ^ 
Tioc, (where r/oc is around lO^^^s at usual tempera- 
tures and densities,) we may increase the range of va- 
lidity of the hydrodynamic electromagnetic theory, from 

^ l/rp to w <C I/tjocj by taking the energy as a local 
function also of P, P, and derive the equation of mo- 
tion for P. There are quite a number of systems with 
a large rp: All electro-rheological fluids have especially 
large rp, of the order of 10~''s, but other complex fluids 
with large molecules and a permanent molecular dipole 
moment (such as nematic liquid crystals) also have a 
slow polarization. (Even the comparatively small wa- 
ter molecule with its permanent dipole moment has a rp 
around 10~^s, just slightly too fast.) The hydrodynamic 
theory of dispersion presented below is unqualifiedly valid 
for these systems, (though it should usually be enough 
to neglect P in the equation of motion, or equivalently, 
exclude P as an additional variable.) We shall refer to 
this scenario as the hydrodynamic dispersion. 

Interestingly, essentially the same set of equations also 
accounts for a system in the ballistic regime, for quickly 
oscillating electric fields and polarizations, l/wp ^ Tioc 
— if we confine our curiosity to questions such as what 
is the averaged force that a high frequency external field 
exerts on a volume element. (Note that the low reso- 
lution is quite sufficient for resolving the hydrodynamic 
responses to a high frequency field.) This scenario we 
shall refer to as ballistic dispersion. Clearly, we need to 
understand why the same set of equations also works for 
the ballistic dispersion, and what the restrictions arc. 

First, a coarse grained, hydrodynamic- type description 
is at all possible because the field variables P, D and 
B vary (due to the largeness of the light velocity) on 
macroscopic, hydrodynamic length scales. Second, most 
of the general principles used as input to consider hydro- 
dynamic dispersion are also valid here. Especially, the 
total energy and momentum remain conserved. The one 
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exception is local equilibrium, the lack of which intro- 
duces some caveats with respect to temporal resolution 
and to dissipation. More specifically, they are: 

• In the ballistic regime, the variables of the theory 

divides into two types, fast and slow. The field 
variables are fast, the rest is slow. They are the 
densities of mass p, entropy s, total energy w, and 
total momentum Therefore, the equations 

of motion of the field variables are highly accu- 
rate, resolving temporal increments much less than 
l/wp, while the actual hydrodynamic equations are 
much coarser, with a resolution low compared to 
Tioc- And because every differential equation, con- 
sistently, has a unique resolution, all field terms ap- 
pearing in the slow, hydrodynamic equations need 
to be appropriately averaged. 

• The generalized hydrodynamic theory presented in 

the next section is shown below to be clearly valid 
for ballistic dispersion in the transparent region, 
where electromagnetic dissipation is negligible, and 
the time rp with which P looses energy long. But 
as argued below, it should remain valid even if field 
dissipation is strong. 

Taking electromagnetic dissipation into account, the 
total, conserved energy divides into three parts, 

U = f/"** + U""^ + C/™^ (2) 

The first is the thermodynamic energy in the absence of 
an external field; the second is the additional energy in 
the presence of a field; and the third is the rest, the en- 
ergy of all microscopic variables not given in the first two 
explicitly, U™'^'^{x\.x^ . . .). The variables .t,; are defined 
such that they vanish in local equilibrium, so they are 
irrelevant for the consideration of hydrodynamic disper- 
sion. In the ballistic regime, is finite and serves 
as a transit hall: External energy is being fed continu- 
ally into U^™, the electromagnetic dissipation excite some 
microscopic degrees of freedom Xi, and convert C/''™ into 
[ymic — which after the comparably long time of Tioc 
becomes heat, C/*"" jTds. The rate at 
which U^"^ is lost is approximately C/''™ « —U^'^/tp, 
the average time this energy stays in the transit hall is 
Tioc, SO U"''" « {tiocItp)^'^. The right side translates 
into {tioc/t)U^'^ for the overdamped oscillation, and into 
^{ujpTiocT)U'^™ for the resonating one. In the first case, 
we always have [7™'^ ^ [/°™, and in the second we mostly 
do, rendering the transit hall usually large. 

Including nonhydrodynamic variables such as P leads 
to contributions ^ P and dU/dF in the energy and mo- 
mentum flux, see next section. If [7™"= is nonzero, we 
would expect similar terms ~ Xi and dU/dxi. These 
we may neglect in the transparent region of vanishing 
dissipation, defined as the frequency regime where the 
stringent condition UpTiocT ^ 1 holds, or equivalently 



jjmic ^ [7°™^ SO and its contributions may be ne- 

glected. Outside these regimes, circumstances arc not 
as certain and in need of a clarifying, more microscopic 
approach such as the Boltzmann theory, to confirm the 
considerations in the next paragraph. 

There are reasons why we may quite generally neglect 
l/mic. -yyjiiig terms such as P^ and P^, of macroscopic 
extent, coherently add up over many periods to yield di- 
verse, slowly varying contributions ~ (P^) and (P^) in 
the momentum and energy fiux, and thereby directly al- 
ter the slow, hydrodynamic variables, the quantities Xi 
are random and of microscopic spatial scales. They fur- 
ther dissipate and degrade, to eventually turn into heat. 
So on a coarse, hydrodynamic time scale, we may sim- 
ply lump (f/'"''=) into heat /Tds, and (^J/^'^) into the 
heat production R. Then, clearly, is neglected as 

an independent entity. 

On a more fundamental level, the very criterion by 
which we have singled out P and P from the lot of mi- 
croscopic degrees of freedom is because they are quali- 
tatively different: Given a certain energy content in the 
field D and B, there is a back and forth of energy flow 
between D, B, P and P within each period; while the 
field energy that leaks into the other microscopic degrees 
of freedom is usually lost. In fact, for an overdamped res- 
onance, it is (as mentioned) appropriate to exclude P as 
an explicit variable, and consider it as one of the many 
ordinary microscopic degrees of freedom, as the energy 
leaked into P is lost to the field. On the other hand, if a 
system involves more variables in the tidal-like transfer 
of field energy, the present theory needs to be general- 
ized to also include them — one example comes readily 
to mind: an independent magnetization. 

Note that being a function also of the slow variables, 
f/*^™ is not conserved by itself, and the permeability e will 
in general contain an imaginary part to account for this 
fact, even without any dissipation (or electric charge). 



B. The Equations of Motion 

The complete hydrodynamic theory of dispersion con- 
sists of a closed set of partial differential equations that 
governs the dynamics of the medium and the electro- 
magnetic field. The structure of the equations is de- 
termined by general principles: the Maxwell equations, 
the Lorentz-Galilean transformation, the thermodynamic 
theory and the relevant conservation laws. 

Combining the two macroscopic energy densities 

we take it as a function of the entropy density s, mass 
density p, the electric and magnetic field D and B, the 
electric polarization P, its canonical conjugate a (that 
will turn out to be essentially ~ P), and the thermody- 
namic momentum density g, 
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^jjM^c ^ _^ ^^^p + vdg 
+H dB + h dP 



E dD 

bda, 



(4) 



where the thermodynamic momentum density g is re- 
lated to the total momentum density 



9*°* = pv 



through 10 



(E X H)/c 



D X B/c. 



(5) 



(6) 



As soon as the energy function U^^'^ is known, the 
temperature T, chemical potential /x, velocity v, field 
strengths E and H are also determined. (In accor- 
dance with iQ], the polarization defined here is a rest 
frame quantity, P = Dq — Eq.) 

Isotropy of space results in the identity 

ExD + HxB + hxP + bxa + ?;xg = 0. (7) 

The Maxwell equations 



VB 
VD 



0, 



B = -cV X E, 
t) = cV xH- p^v 



(8) 



account for the motion of D and B. Here, the dot indi- 
cates partial temporal derivative d/dt and denotes the 
macroscopic charge density. The variables p, U, are 
conserved, their equations of motion take the form 



p + V- (pv) = 0, 
f7-HV-Q = 0, 



(9) 
(10) 

(11) 



where Q is the total energy flux, and (n.y — 11^) the 
symmetric total momentum flux, or total stress tensor. 
The entropy is not conserved, and has a positive source 
R, 



s + V ■{sv-{^) = R/T, R>0, 



(12) 



The dissipative part of entropy flux describes espe- 
cially thermoconduction, while Hf^ accounts primarily for 
viscosity-related phenomena. 

The macroscopic variables P and a are governed by 
equations that are essentially of the Hamiltonian type, 

P = dU/da = b, a = -dU/dP = -h, 

with some supplementary terms needed to ensure the 
proper transformation behavior, and to account for dissi- 
pation. First, the temporal derivative is replaced by the 
Galilean invariant operator that takes into account the 
effect of the local movement of the medium, 



Dt = dt + ivV)-nx, 



(13) 



where fl = x v. Second, a dissipative force h-^ is 
introduced in the equation for a to account for electro- 
magnetic dissipation that (in the linear case) is usually 



taken care of by an imaginary term in the electric permit- 
tivity e. Third, the polarization is changed if the medium 
undergoes volume dilatation, as a term — P(Vi') appears 
in the equation of motion for P, 



AP = b-P(Vt;), Aa=-h-h^ 



(14) 



The dissipative terms are determined with the familiar 
method of irreversible thermodynamics: We first identify 
the entropy production as 



R + J7""' = f ^ • (VT) + • b + n 



ij '^V ' 



(15) 



where Vij — ^(VjUj -I- \7jVi). Then take the fluxes as 
proportional to the thermodynamic forces, 



/f 1 = I A,, ) X ( V,r ) (16) 



^tjl ^IJ ^Hj 



(Appropriate Onsager reciprocity relations are implied.) 
The energy flux is 

Q, = {Ts + tig + Vkgk)v^ - TfP - v,Uf^ + c(E x H); 

(17) 



-1- v,{h - P) + ^[v X {hxP + h X a) 



(The last two terms were erroneously omitted from [Q.) 
This expression may be rewritten as 



= c(E" X H'J), - //^T + C/^^^-^z;, 
+ (n,, - nP^)vj - Vkgl°'v,, 



(18) 



to see that the velocity-dependent terms do come from 
an Lorentz-Galilean transformation, discussed eg in 
(E° = E + t) X B/c and HO = H - v X D/c arc the 
restframe flclds.) 

The stress tensor is symmetric and given as 

Ily = ^[vigj - EiD-j - HiB-j + (i ^ j)] + {Ts + pp 
-h£/-t;-fH-B + E- D-|-h-P - U^^'"')5.,j. (19) 

Frequently, there are many different resonance frequen- 
cies of the polarization, not just the single one, given here 
as LOp. This fact can be accounted for by introducing as 
many "subpolarizations" , 



D E = P = ^ Pa, 



(20) 



chosen such that the two squared order terms of the en- 
ergy are diagonal. 



rrcm _ 



(21) 



Close to one resonance a, if it is well separated, as all 
the other subpolarizations are not excited, we may sim- 
ply substitute Pq, for P. 
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C. Some Explicit Expressions 

Now, the above equations are rendered more explicit 
by an expansion of the energy function in the vector- 
variables D, B, P, a and v to third order, as this is suf- 
ficient for a comparison to the linear results by Brillouin 
and Pitaevskii's. For a magnetically inactive medium (ie 
taking the static magnetic permeability as 1), such an 
expansion yields 



U 



Ma 



IB- 



Id'- 



- • (P X a) + \pv- + 



D-P + ipVx 

^pv- + 0\ (22) 



where O" denotes terms of n-th or higher order in the 
vector variables (D, B, P, a, d). The energy density in 
the absence of electromagnetic fields is ?7™'**(s,p); the 
coeSicient x is related to the conventional static dielec- 
tric susceptibility x' = Pl^ by x'^^ = X^^ ~ 1; '^p is 
the dielectric resonance frequency; ^ is connected to the 
magnetic cyclotron- frequency ws; all these parameters 
are in principle functions of p and s. 

Obtaining the differential form from Eqs(^, ^) 



that may be used to eliminate a in the above formulas, 
and write them instead with P. Especially, the mag- 
netization in the rest frame and to lowest order is then 
(f/x'^p)P X P. The nonmagnetic magnetization is now 
seen to result from rotations of the polarization. 

For a qualitative estimate of the coefficient envi- 
sion electrons revolving around their ion centers |^ . As- 
suming the rotations to occur in phase, the magnetiza- 
tion associated with it is ((7e'^e/2'Tie)L, where and 
are the charge and mass of the electrons, while rie de- 
notes their density.The angular moment of the electrons, 
L = TOefe X i"e (with the radius of the circular motion) 
can also be written as L = {ra^j q-^iJ^V x P, because the 
polarization P is qeneVe. The attendant magnetization 
is M = (l/2(7ene)P x P. So the coefficient ^ is 



(31) 



1, and LUp can 



Particularly for an electron plasma, x 
be considered as the plasma frequency {q1nf,/ra^Y/- . 
Eq(^) reduces to ^ = qe/^nie — —ujb/'^B, with ujb = 
—Bqe/nie the plasma cyclotron frequency. 



^^f^Mac _ ^ . = + pdp + E dU + H dB 



III. MONOCHROMATIC APPROXIMATIONS 



h • dP + h ■ da — g ■ dv, 



(23) 



we can derive the thermodynamically conjugate vari- 
ables, 



T = 



ds 2 ds 
OS 



2 ds 



(24) 



dp 



2 dp 2 dp 







(25) 


E 


= D-P-t;xM + C'3, 


(26) 


H 


= B - ^(P x a) + it; X P + ©3 

c 


(27) 


h 


1 1 o 
= -P - D - ^a X B + -B X t) + 0^ 


(28) 




X c 


b 


= Xo;pa-^B X P + C^ 


(29) 



Note that the magnetization M = B — H as given 
in Eq(p7|) is a term of order . So the difference in 
the polarization v x M, between the rest frame quantity 
Do — Eq and the laboratory quantity D — E, is of order 
. Within the accuracy of the above equations, it is 
therefore ignored. 

Inserting ( p9[ ) in (|l^), we obtain the expression 

a = ^ [AP + P(V-i;) -f CB x P] -f ©3 (g^^ 
X'^p 



With the help of a closed L-C circuit, Pitaevskii ob- 
tained a number of important results on the effects of 
a high-frequency field in a medium Q. Because of the 
special setup, the results are subject to certain restric- 
tions. In order to compare our theory with his work, the 
same limits will be taken in our theory. Therefore, we 
shall consider a transparent medium exposed to a strictly 
monochromatic electric field: 



E 



1 



-I- C.C., 



£ = 0, 



(32) 



where £ is the constant amplitude and lo the frequency. 
From now on, we shall always assume that the medium 
is at rest, so any velocity-dependent terms will be dis- 
carded. Because only the electric properties are of inter- 
est, we also omit the quickly oscillating part of magnetic 
field in the medium, as in though a strong, constant 
magnetic field is allowed to be present. If the material 
coefficients x^ ^pi ^ constant with respect to time, 
the induction D and polarization P will also take the 
monochromatic form 



D = ipe"''^' 
2 



c.c, P = -Ve' 



+ c.c. 



(33) 



However, if Xi '^pi C ^^'^ B are allowed to vary - slowly 
- via their dependence on the density or temperature, 
the fields D and P will become quasi-monochromatic. 
In this case, we may still hold w to be strictly constant, 
while allowing the amplitudes V, V to chang slowly with 
time. The qu asi- m onochromatic situation will be studied 
in subsection HIE. 
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In what follows, the dynamic equations given in the 
previous section will be investigated, under the precondi- 
tions mentioned above. We will show in detail the deriva- 
tions of four formulas, all well-known in the literature [Q. 



A. The permittivity 

The frequency-dependent permittivity Sij is calculated 
from the equation of motion for P. The expression is 
given by inserting (|3^) and (|2|) to (|l^) , taking the coef- 
ficients X, ^ and the magnetic field B as constants. Ne- 
glecting the velocity-dependent terms, we have 



_ P £B X P P 

D 2^ ^ = 0. 



(34) 



If the fields E, P assume the monochromatic form of 
Eqs(p3,R3h, the above equation becomes 



1-- + ^ )P + 2*^B xT' + f^O. (35) 
X X^^p J X^p 



Solving it for V, we obtain 

P=(£i-l)£: + e2(B£:)B-Kz£3(£:xB), (36) 

with 

XUJpiuj^ ~ ujp + X^p) 



£1 = 1 



(37) 



£2 



and 



where 



S3 



XUJpUJUJB 



CUB = -2^B. 



(38) 



(39) 



(40) 



Using the fact that Pi = {£ik — 5ik)Sk, we observe from 
(Ba) that the permittivity is 



£iSik + SiBiBk + iss^ikiBi 



(41) 



where eiki is the total antisymmetric tensor, £123 = 1. In 
the low frequency limit w ^ 0, we have according to (37) 
El = 1/(1 — x)- Note that the imaginary term in (41) 
is not connected to dissipation. It is a purely reactive 
term. This can best be see from its invariance under the 
time-reversal operation: lu —f —uj, Bi — )■ —Bi. 



B. The energy density 

Eliminating the quantity a in the energy function (|22 
with the help of (|3C|), we get, for a medium at rest and 
including terms of third order in the field 



B^ 



p2 

DP + — 



2X 2x^2 



(42) 



Now consider the monochromatic case (p2,p3) and apply 
a time-average procedure denoted as (^-77 the energy 
density is then given as 



52 



1 



([/) = + — + -£k£l + 



1 

4x 



1 



- X VkVl 



or 



2 4 



- - 

4 \x 



X^l 



1 ) (£fem — 5km){£*kn ~ Skn)£rn£n^ 

(43) 



where * denotes complex conjugation. Because in our 
work the permittivity is given by ( p| ) and (|37|-p9|), we 
can verify by direct computations that the equation 



duj 



^mn ^" [1 

X 



X {el 



H n ~ X] [Skm — 5km) 

^p J 

Skn) (44) 



holds for this form of permittivity. So the time-averaged 
energy density could be expressed as 



n2 

{U) = U"^"' + ^ + 



1 duJEmn 

4 duj 



F F* 



(45) 



which is the Brillouin's expression for the time-averaged 
energy density of the electric field ||] . 



C. Pitaevskii magnetization 

The equation ( p7| ) shows that a magnetization could 
be induced dynamically in an electrically polarizable 
medium, although the static magnetic permeability is 1. 
Inserting (30) in ( p7| ) we obtain, for a medium at rest, 
the magnetization 



M 



X^p 



■P X P + O^ 



(46) 



In the monochromatic approximation, the time-average 
of this magnetization is the same as that obtained by 
Pitaevskii H. Indeed, inserting (^, p3|) in (E6|) wc have 
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iijkiejm - Sjm){el„ - 5kn)Sra£l- (47) 



^ ■' op 



With the help of the expression of the permittivity (41) 
Ipy) , one can show the vaHdity of the equation 



2 dB, 



So (|3) possesses the form given by Pitaevskii 



(M) = 



4 9B 



It-* p 



(48) 



(49) 



Axi^lH^ dp i ' 

- ^-{Z.{D] + Epi + c.c) - i((Af),i/, + (A/),i/0- (52) 

Because the difference between B and H (ie the 
Pitaevskii magnetization ([49|) ) is of second order in f , 
we can write the formula (pq) for V to the same accu- 
racy as 



D. The stress tensor 



7'=(£i-l)f + £2(Hf)H + i£3(f xH), (53) 



Inserting the expressions (gg), (|3|), (|2|) for [/, T, ^ 
into the stress (|l9|) and eliminating the fields B, a, h 
with the help of (p7|), (28), (|30|), the stress becomes 



n, 



-^H P X P + — p—- + -^^p^-J^ 

' 2 V dpx) ^^^p dp 



^H(P X P) p 



AX 



-[E,D,+KUH, + {i^j)]+0\ 



(50) 



with 



£1 = 1 



(^2 - ^2 + Xt^p)^ - ^^^^^ff ' 



£2 



and 
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£3 



XUJpUJl^H 



(54) 



[(u;2 - w2 + ^^^2)2 „ tj2w2^](tj2 - ^ X^p)-^"^ 

(55) 



[(tj2 - Cj2 + xtj2)2 _ uj'^uj%]H' 



(56) 



where 



where po is pressure of the medium in the case without 
electromagnetic fields, 



poip,s) = -U'^^' + p 



dp 



ds 



(51) 



(57) 



Now inserting in the stress (|5^) the expression ([5^), it 
becomes 



Here, in order to avoid unnecessarily complicated formu- 
las, we also neglect the entropy dependence of the pa- 
rameters Xi ^pi C in (Fo|)- When the two electric fields 
take the monochromatic form (^ ^3|), we obtain after 
the time-averaging procedure. 



(n„> = { 



tt2 1 



uj^ ~ujp+ X^p 
4x'^p 



■p-p* 



(n,,> = { 



i?2 1 / de 



H,Hj - hs.V* + £,V* + C.C 



-{{M),H, + {M),Hi), 



(58) 



2xwi 



where we have used the following equations valid for the 
£i, £2, £3 given by (|4[j~ 
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dei _ujI^^ 2uju;hH ^ 

^ - - 1) + ^3^ ] - — ^^3(£l - 1), 



2ujluh 



+ 



[2s2{e I - 1) + siH' - ei] 

,7 TT 

esiei - 1), 



X^f^lH 



(59) 



(60) 



"g^ - ^ .,2. .2 ^3(^1 ^ 1) " .,2. .2u i(.'^l - 1) + J. 



2 2 '^3('^1 ~ 1) ^ 2 2 1/1 



(61) 



1^3 



[f *(^^ X H)j + £»(H X + H,(r X £).j\ 



+ {i^3)]. (68) 

If the following identity is noted 

£*(£: X H)j + £,(H X £*)j + i/i(£:* x £)j 

= {£*xE)m^j, (69) 

we finally obtain 

£T Em ( dSlm £lm ~^ st 



4 V dp 



01' 



dei 



2luujhH 
X^p 



e3(ei - 1), 



+ 



2ujujh 

x4^ 



£3(21 - 1), 



(62) 



(63) 



-^[sAS* + e2(H£*)(£,i/j + £j£fO + c.c] - H,H, 



+ 1 



dei 
'dip 



EE 



dm 

{E* X E)ii 



(H£:)(Hf*) 



(70) 



This agrees with the Pitaevskii's stress tensor of a vari- 
able electric field in a liquid located in a strong magnetic 
field il. 



^£3 



2u?_ 

X^t 



£^{£1 - 1) - 



X4H 



dei 2uH . ^. 
= 2"'^3(ei - 1), 



dujH x^p 



de2 2uj 

aJ7 = -;;7T^^3(£i-i), 



H X^pH 



dujH X^jH 



(64) 



(65) 



(66) 



(67) 



Using ( [49D and the fact that Vi — EimEm, we can also 
write the tensor (58) into the form 



El Em ^ ^ ds Irn 

dp 



— £lrn — HiHj 



-HiHj 



\{£AE; + S2{B.E*)E,H, + e2{liE)E:H, 



^^^EE* + ^(Hf )(H£:*) +z^(r X £)H 



E. Time-dependent permittivity 

In the previous subsections, we assumed that the pa- 
rameters X, ojp, £, and the magnetic field B are time- 
independent. Co nsequ ently, the permittivity discussed 
in the subsection III A is static. Now, we will abandon 
the restriction and allow x, i^p, B to change slowly with 
time. This case is naturally accounted for by (p^, ( [2^ ) 
and (|30|), the equation of motion for P, 



d_ 

dt 



X^l 



2^B 

X^l 



1 

1 

X 



X P - 



d gB 

dt xLoI 



X P 



P-E = 0. 



(71) 



Here, we again neglected the dissipation and con- 
sidered a stationary medium: v{r,t) = 0. Comparing 
the equation (|7l| ) with (^J), we can see that the tempo- 
ral variations of Xi'^pi'CjB give rise to additional terms, 

which result in a dynamic correction efj" to the static 
dielectric permittivity obtained in subsection III. A. In 
other words, the relationship between the amplitude of 
polarization V and that of electric field E is no longer 
given by (pq), but by 
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'Pi — (Sij 



dyne 



(72) 



where is given by ( ^ ) and ( ^7[]39| ). e^J" may be calcu- 
lated by retaining a monochromatic electric field in (71): 
£ = 0. Yet, because e is now time-dependent, the ampli- 
tude of polarization V will change with time, Vi — e 
And Eq([7l|) becomes. 



IV. DISCUSSIONS 



1 ^ \ S.^ ^ 

1 - - + — ^ + 2i-^B X V 

X X^p ) X^p 



-£ + ^ \ — -e£ X e£ 



— lU) 



dt \x^^p X^l 
^ ^ fl-^ ) X (£ + l)£ = 0. 



dt xu^l 



dt xL^p 



where e is matrix notation of the permittivity 
ing for V, we obtain 

2iuj . . f d I 



(73) 
Solv- 



2 £mj + «CJ 



dt x^^l 



2^ 



iei,-di,) . (74) 



x^^p " ydtx^l, 

Together with the static permittivity derived before, this 
dynamic correction provides the expression for the full 
time-dependent permittivity. 

To compare with we decompose the dynamic cor- 
rection d?^) into a hermitian and an antihermitian part: 



dyn _ / dyn , ■ // dyn 



(75) 

Both matrices 'e'^y" and "e'^y" may have complex ele- 
ments, but must be hermitian. In accordance to (|75|), we 
may also call them the real and imaginary parts of e 
Particularly, the imaginary part is 



dyn 



"j-dyn _ / dyn* 



-dyns. 



-V 2 

which can be also written as 



d_ 

dt 



x^l 



(76) 

here the fact e*j = Eji is used. From the expression ( ^ ) 
and (|37|-^9|) for e, one can show that the equation 

(77) 

is valid. So the imaginary part ( fz^ ) of the dynamic con- 
tribution to the permittivity is 



1 

//^dyn _ ^ ^ ^tj 



''^ 2dujdt' 
This formula was first obtained in 



(78) 



Because both the dispersion and nonlinearity are ac- 
counted for, the hydrodynamic theory of dispersion 
sketched in section II is a fairly complete theory for 
the dynamics of a fluid interacting with varying fields. 
The theory is derived by generalizing the hydrodynamic 
approach, but the result is consistent with the work of 
Pitaevskii, who starts from rather different physics. 

Though not shown here, the present theory reduces to 
the hydrodynamic one |^ , in the low frequency limit lu — > 
0. It is also in agreement with the Barash-Karpman's 
stress tensor derived for quasi-monochromatic field (ie 
including the lowest-order effects of temporal variation 
of the field amplitude £) ^ . All these features support 
the statement that the basic equations shown in section 
II are correctly formulated, particularly the fundamental 
differential relation (Q). 

In our theory for dispersive media, the Pitaevskii's 
magnetization appears as a consequence of circular mo- 
tions of the polarization. In contrast to the conventional 
magnetization of atomic origin, the Pitaevskii's magne- 
tization is macroscopic. Because the circular motion of 
P is usually accompanied by a rotating electric field, the 
Pitaevskii's magnetization is less suitably generated by 
linearly polarized electromagnetic fields, as suggested in 
It is remarkable that the explanation to this phe- 
nomenon was given by the first nonlinear term in the 
expansion of the energy, Eq(|2^). We expect it to be an 
important nonlinear effect, and to play a significant role 
in the nonlinear interaction between matter and intense 
laser lights. 
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